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ABSTRACT 

We consider the inward propagation of warping and eccentric disturbances in discs 
around black holes under a wide variety of conditions. In our calculations we use 
secular theories of warped and eccentric discs and assume the deformations to be 
stationary and propagating in a disc model similar to regions (a) and (b) of Shakura 
and Sunyaev discs. We find that the propagation of deformations to the innermost 
regions of the disc is facilitated for low viscous damping and high accretion rate. 
We relate our results to the possible excitation of trapped inertial modes, and to 
the observations of high-frequency quasi-periodic oscillations (QPOs) in black hole 
systems in the very high spectral state. 

Key words: accretion, accretion discs — black hole physics — hydrodynamics — 
waves — X-rays: binaries 



1 INTRODUCTION 
1.1 Background 

In the classical theory of accretion discs the gas is assumed 
to rotate around a massive central object following copla- 
nar, circular and approximately Keplerian orbits. Although 
this is the simplest solution for the motion of fluid elements 
placed in a Newtonian potential well, external forces can 
give rise to non-planar discs, composed of non-circular rings. 
In fact, the general solution allows for discs to be twisted or 
tilted, with their orbits presenting a smoothly varying eccen- 
tricity and/or inclination. The litera ture offers strong ob- 
servational and theoretical evidence (IBardeen fc Pettersonl 



1 19751 ; IOgilvidl2000l . |200ll : iFragile et al"1l2007l , and references 
therein) to believe that under various conditions discs can 
become either warped or eccentric. In particular, the pre- 
cession of deformed discs has bee n used to explain several 
phenomena in X-ray binaries (e.g. iGerend fc Boynton|[T976l ; 
IWhitehurstlll988l ; IStella fc Vietri|[l998l ; lLaj|l999h ~ 

Discs around black holes are warped if the spin axis 
of the central object and the angular momentum of the 
accreting matter are misaligned. The fluid elements in 
tilted orbits will be subject to Lense-Thirring (gravitomag- 
netic) precession which ten ds to twist and warp the disc 
|Bardeen fc PettersorJll975l l. The way a tilted ring commu- 
nicates its inclination with respect to the equator of the 
black hole to its neighbouring rings depends on a balance be- 
tween hydrodynamic and gravitomagnetic torques. Roughly 
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speaking the warp propagation is diffu sive if a > H/r, and 
is wav e-like if a < H/r, where a is the lShakura fc Sunvaevl 
(|l97&t ) viscosity parameter and H/r is the angular semi- 
thickness of the disc. In a low-viscosity disc around a Kerr 
black hole, t he warp has an osc i llator y radial structure, 
as shown by llvanov fe Illarionovl (|l997l ) and discussed by 
iLubow et al. I (|2002h . In several X-ray binaries, jets (which 
are thought to be aligned with the rotation axis of the black 
hole) are observed to be misaligned with respect to the b i- 
nary rotation axis (e.g. GRO J1655-40. [Maccaronel 120021 ) ■ 
This is strong evidence for a warped disc around the com - 
pact object in the binary (|Martin. Tout fc Pringld 120081 ). 
Warping migh t also be induced by radiation pressure forces 
|Prmglelll996l ). Despite being more likely for d iscs around 
neutron star primaries (|Ogilvie fc Dubual200ll '). radiation- 
driven warping might also occur around black holes pro- 
vided the disc exte nds to large enou gh radius, as in the case 
of GRS 1915+105 (]Rau et alj|200l) . 

Eccentric accretion discs are believed to exist in cata- 
clysm ic binaries of mass ratio q < 0.3 (e.g. iPatterson et al.l 
l2005h . Eccentricity can result from an instability involv- 
ing the orbiting gas and the tidal potential of the com- 
panion star. Discs may become eccentric if they are large 
enough to extend to the 3:1 resonance, the radius in the 
disc where i ts angular vel ocity equals three times that of 
the binary jLubowl [i991al lbT). The strongest observational 
evidence for eccentric discs is the detection, in the light 
curves of accreting binary systems, of long-period modu- 
lations known as superhumps, which can be explained by 
the action of tida l stresses on a precessing eccentric disc 
l|Whitehurstlll988T ) . Although this phenomenon is usually as- 
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sociated with cataclysmic variable stars, superhumps have 
been d etected in an increasing num b er of black hole bi- 
naries dOTlonoghue fc Charles) Il996l: lHaswell et alJ l200ll : 



lllemura et al.ll2002l ; iNeil et al.H2007l ; Izurita et al.ll2008l ). In 



fact, black hole binaries are likely to have mass ratios q < 0.3 
and therefore to have ecc entric discs during; a t least some 
stages of their outbursts. lHaswell et alJ (|200ll ) noted that 
the mechanism responsible for superhump luminosity varia- 
tions differs in X-ray binaries and cataclysmic variables, but 
the underlying dynamics is the same. 

Our goal is to describe the propagation of global eccen- 
tric and warped disturbances, under a variety of conditions, 
from the outer parts of the disc where they originate to the 
inner parts. In the outer parts of the disc, the motion is 
Keplerian and stationary or slowly precessing global defor- 
mations are supported. Since the Lense-Thirring precession 
frequency increases with decreasing r, tilted orbits at differ- 
ent radius will tend to precess at different rates, which tends 
to twist the disc. Hydrodynamic stresses counteract this ef- 
fect, resulting in wavelike warp propagation from the outer 
disc to the inner region. Similarly, for eccentric discs a tran- 
sition has to be made between the Keplerian region where 
eccentric instabilities are driven, and the inner region where 
general relativistic effects dominate the precession of ellipti- 
cal orbits. It is this connection between different regions of 
the disc, with different characteristics, that we attempt to 
describe. 



1.2 Deformations as global modes 

Warped or eccentric disturbances of small amplitude can be 
thought of as linear perturbations of a circular and coplanar 
disc. In the simplest case of a vertically isothermal disc with 
unit ratio of specific heats, wave modes with a dependence 
in time and azimuth of the form exp (im(f> — iuit ), and local 
radial wavenumber k obey (|Okazaki et al.lll987f ) 



|2 (^-K 2 )^ 2 



nVLi 



(1) 



where ui = u — mQ is the Doppler-shifted wave frequency, 
integers m and n ^ are the azimuthal and vertical mode 
numbers, respectively, and c s — Q. Z H is the constant sound 
speed in the disc, with H being the local vertical scaleheight. 
The symbols f2, n and Q z represent the angular velocity, 
epicyclic frequency and vertical frequency, respectively. Rel- 
ativistic effects are included in the problem by using rel- 
ativistic expressions for these characteristic frequencies. In 
the Kerr m etric, valid for a rotating black hole, they read 
|Katolll990l ) 
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where a is the dimensionless spin parameter of the black hole 
(— 1 < a < 1), r is the cylindrical radius in units of GM/c 2 , 
and the frequencies are in units of c 3 /GM, where c is the 
speed of light, G the gravitational constant and M the mass 
of the black hole; in these units velocities are in units of c. 
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Figure 1. Local wavelength of (a) (u>,m,n) = (0,1,1) and (b) 
(0,1,0) modes (stationary warp and eccentricity, respectively) in 
units of H for several values of the black hole spin. The inner 
edge of the disc is taken to be the marginally stable orbit, where 
k = 0. 



The dispersion relation is local, i.e., it is different at each 
radius since the characteristic frequencies depend on r, and 
is valid only where the wavelength A = 27r/fc <C r. 

Viewed as linear perturbations of a standard disc, warp 
and eccentricity correspond to low-frequ ency waves w i th az - 
imuthal number m = 1. As showed by iKatol (|l983l . Il989h 
these modes are global, and have relatively long wavelengths 
when compared to the other modes described by the local 
dispersion relation, being the most likely to exist in a tur- 
bulent disc. More specifically, a (m,n) — (1,1) mode can 
be identified with a warp as its vertical displacement is in- 
dependent of the vertical coordinate z and varies with <f> at 
each radius, giving rise to a tilting (movement with respect 
to the disc plane). A (m,n) = (1,0) mode corresponds to 
eccentric disc since its radial displacement is independent 
of z and varies with azimuthal angle, originating an elliptic 
orbit (movement within the disc plane). The m = 1 modes 
are global in the sense that only they can have a wavelength 
much larger than the semi-thickness of the disc over a wide 
range of radius. This can be confirmed in Fig. [T] where we 
show the variation with r of the local wavelength in units of 
H, for several values of a, for bot h the warp and eccentricity. 

In I Ferreira fc Qgilviel (|2008l ) we worked with a simpli- 
fied set of hydrodynamic equations in order to deal with 
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the nonlinear couplings of different wave modes. Within this 
framework we obtained an equation for t he warp tilt W at 
each radius (similar to the one given by IPapaloizou fe Linl 
1 19951 ) assuming the warp to be a zero-frequency mode prop- 
agating in a strictly isothermal, relativistic disc. In addition 
to being valid only where c s is constant and where viscous 
effects are negligible, the equation derived does not hold 
for large radius, where the wavelength of the warp becomes 
comparable to r, since some terms were omitted in our lo- 
cal approximation. A similar method was used to derive an 
equation for the eccentricity at each radius E valid only for 
an inviscid and strictly isothermal disc. 

Here we consider the more general description of 
a stationary , wave - like warp and eccentricit y give n by 
iLubow et all (|2002f ) and IGoodchild fc Ogilviel (|2006l ), re- 
spectively. These th e ories are les s general than the ones for- 
mulated by I Ogilviel (|2000l . l200ll ) since they consider defor- 
mations to have small amplitude, and not all viscous or tur- 
bulent effects are included. Also, the IGoodchild fc Ogilviel 
|2006l ) calculation uses a highly s implified, 2D disc m odel. 
However, the secular theo ries of ILubow et al.l (|2002T l and 
IGoodchild fc Ogilviel \ 20061 ) are appropriate for thin discs of 
any given structure, and allow for viscous (turbulent) damp- 
ing of deformations. 



1.3 Black hole states and high-frequency QPOs 

Observations indicate that black hole systems can be found 
in different spectral states, going from the quiescent state to 
low, intermediate, hig h and very high states as the ma ss ac- 
cretion rate increases |Remillard fc McClintockll2006al ). The 
very high state is the one where high-frequency QPOs, i.e., 
modulations seen in the spectra of black hole systems that 
are not strictly periodic and have frequencies > 100 Hz, are 
almost exclusively detected 0- This is the dominating state 
close to the Eddington limit. 

The accretion flow has different characteristics in dif- 
ferent black hole states. If the disc properties are such that 
the warp and eccentricity can propagate to the inner re- 
gion, the interaction between relativistic effects and global 
deformations can give rise to interesting phenomena. In par- 
ticular, the warp and eccentricity can have a fundamen- 
tal role in exciting trapped i nertial modes, whi ch may ex- 
plain high-frequen cy QPOs ([Nowak et al] 1997j~) . as previ- 
ously suggested bv lKatol (|2004l . 120081 ) and lFerreira fc Ogilviel 
|2008l ). 

In this paper we show that the inward propagation of 
warp and eccentricity is facilitated for high accretion rate. 
For fixed viscosity, global deformations reach the inner re- 
gion of the disc with a modest amplitude when the accretion 
rate is high, and can take part in the excitation mechanism 
for inertial modes trapped in this region. If high-frequency 
QPOs can be identified with these trapped inertial waves, 
we can argue that they are predominantly detected in the 
very high state, where the accretion rate is close to Edding- 
ton, because only in this case global deformations can reach 



the inner region with non-negligible amplitude and excite 
the trapped modes. 



1.4 Plan of the paper 

In section 2, we introduce the disc model in which we study 
the propagation of global deformations. We consider a disc 
with a poly tropic structure in the vertical direction, con- 
stant opacity, and with gas and radiation contributing to 
the total pressure. In section 3 we describe the propagation 
of global disturbances, by solving equations for the warp tilt 
and eccentricity at each radius. To mimic several possible 
conditions in a disc with a fixed viscosity a around a black 
hole of a given mass and spin, we vary the accretion rate and 
the damping of the warp and eccentricity. The results ob- 
tained are discussed in section 4. Conclusions are presented 
in section 5. 



2 DISC MODEL 

2.1 Viscosity prescription 

In the famous IShakura fc Sunvaevl (|l973l ) model, the disc is 
assumed to be geometrically thin and optically thick. The 
gas in the disc is heated by viscous dissipation and cooled by 
radiative diffusion through the vertical boundaries. Angular 
momentum is transferred outwards at a rate which is taken 
to be proportional to the height-integrated total pressure, 
with the efficiency of the process being parametrised by a 
dimensionless quantity a, independently of the nature of 
the stress. The disc is taken to be composed of three parts: 
in the innermost regions (a) and (b), electron scattering is 
the main contribution to the opacity (which is then given by 
the constant, Thomson value, kt), while free-free absorption 
is dominant in the outer region (c). In (a) the pressure is 
determined by the radiation pressure (p tx T 4 ), while gas 
pressure (p oc pT) dominates in (b) and (c). The radiation- 
dominated region exists only if the accretion rate is larger 
than a few percent of the Eddington rate, being important 
in accretion discs around black holes. 

Despite neglecting to treat the vertical structure of 
the disc, not considering other contributions to the cool- 
ing and heating rates of the gas, and despite the sim- 
ple parametrization for the angular momentum transfer, 
the Shakura & Sunyaev a-model is still the standard ref- 
erence for the radial structure of accretion discs, if gas 
pressure is dominant. Region (a) has always been more 
problematic. In fact, the assumption that stress scales 
with the total pressure in th e radiation pressure domi - 
nated regime leads to viscous (Light man fc Eardlevl ll974T ) 
and thermal (|Shakura fc Sunvaevl 19761 ) instabilities. These 
can be avoided by sim ply assuming that the stre sses scale 



1 In the X-ray binary XTEJ1550-564, a HFQPO was detect ed at 
a photon count rate below the VHS (see Ivan der Klisll20ul and 
references therein). 



can be avoided by simply assuming th at tne stre sses scale 
with the gas pressure (|Sakimoto fc Coroniti|[l98ll ). or with 
a combination of both gas and radiatio n pressure (e.g. 
iTaam fc Linl 1 1984 iMerloni fc Fabianll2002l ). recent simula- 
tions by Hirose et al. (ApJ, submitted) suggest that turbu- 
lent stresses may indeed be correlated with the total pres- 
sure. However, fluctuations in stress (no w known to be mag- 
netic in nature, iBalbus fc Hawlevlll99ll ) are the ones driving 
pressure fluctuations, contrary to the usual assumption that 
changes in pressure give rise to changes in stress; this seems 
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to avoid thermal instability (which has a faster growth rate 
than viscous instability). With these results in mind, we 
consider the a viscosity prescription in the form, 



H — a 



(5) 



so that the viscosity p and total pressure p correlate, as 
in the original Shakura & Sunyaev model. Here we con- 
sider a to be a function of radius only; in this section 
SI = Qa" = ■J GMjR? is the Keplerian angular velocity, 
and R — rGM/c 2 is the dimensional radius. 



2.2 Vertical structure 

Independently of the viscosity prescription, the radiation- 
pressure dominated region might suffer from another type 
of instability. In the Shakura & Sunyaev model, the dissipa- 
tion rate per unit volume is independent of z. If the dissipa- 
tion rate per unit mass is also vertically constant, then the 
density is independent of z, vanishing a bruptly at the verti- 
cal b oundaries of the disc, z = zq(R) (|Shakura fe Sunyaev! 
1 19731 ). A di sc with these properties is subject to convective 

instability (|Bisnovatvi-Kogan fc Blinnikovlll977r ). 

Simulations of radiation-dominated discs by I Turner! 

(j2004l ) indicate that neither is the dissipation per unit mass 
constant nor is the density independent of z. In fact, the 
computed density profile resembles more that of a po lytropic 
model with index s — 3 (see also lAgol et "al Using 
the polytropic law, p oc p 4/3 , in the hydrostratic equilibrium 
equation 

dp 

dz 



-n 2 z pz (Keplerian disc: 17 z = Q. = \J GM/R S ), (6) 



it is straightforward to get 

3 



P = po(R) 



1 - I — 

zo 



( — 20 < Z < Zo). 



(7) 



This result was used by iBisnovatvi-Kogan fc Blinnikovl 
(1977). who argued that convection in the radiation-pressure 
dominated region would establish a vertically isentropic 
structure {T 3 /p — constant in z). 

The s — 3 polytropic structure has the convenient prop- 
erty of allowing for a ratio of radiation to gas pressure, 



p T _ 4<r^ m m p T 



3cfcs 



(8) 



independent of z. Using Q in the hydrostratic balance equa- 
tion we get p = po(R) [l — (z/zq) 2 ] , where 



Po = 



fl^Zp 

8 



-po = 



Q?zl 



-po- 



This is consistent with the equation of state 



P = Ps + Pr 



-^pT + ^T\ 
/i m m p 3c 



(9) 



(10) 



provided T = Tq(R) [l — (z/zo) 2 ] , and independently of the 
variation of /3 with radius, i.e., the polytropic vertical struc- 
ture with s = 3 can be used to model not only the radiation 
pressure dominated regime but also the gas pressure domi- 
nated region. In fact, a full treatment of the vertical struc- 



ture of a gas pressure dom inated disc in hyd rostatic ((6} and 
radiative equilibrium fe.g. iFrank et alj|2002h : 



^— = fi (i?S7') 2 (energy balance) , 



(11) 



F — F ra d = — (radiative diffusion law), (12) 

3ktP oz 

indicates that the resulting profiles are very similar (ex- 
cept perhaps F(z) for large z) to those obtained for a poly- 
tropic model with s fn 2.7 (Ogilvie 2001, unpublished lecture 
notes). This result, and the simulations in radiation-pressure 
dominated regions, indicate that assuming a polytropic ver- 
tical structure with s — 3 throughout the radial extent of the 
disc is a good approximation, provided the opacity remains 
constant. This is true of the innermost regions ((a) and (b) 
in the Shakura & Sunyaev model) of the disc. Since we are 
most interested in the propagation of global deformations in 
these regions, we take the opacity to be constant throughout 
the disc, i.e., region (c) is neglected in our model. In fact, the 
physics of region (c) is mo re complicated than assumed by 
IShakura fc Sunvaevl (|l973l ). because of the presence of dif- 
ferent sources of opacity and the importance of irradiation 
and partial ionisation. 



2.3 Radial structure 

The assumptions of hydrostatic equilibrium and vertical 
isentropy (or equivalently, s = 3 polytropic structure in 
z), together with the equation of state, are enough to de- 
termine the vertical structure of the disc, i.e., determine 
the variation of p, p and T with z, without the require- 
ment of radiativ e balance (see discussion in section 3.1 of 
lAgol et al.ll200lf ). However, to fully determine the radial 
structure: po(R), po(R), To(R) and zo{R), two more equa- 
tions are needed since the polytropic relation is not valid in 
the radial direction. 

Since matter is being accreted in the disc, there are 
radial drift motions which transport mass and angular mo- 
mentum. In steady thin discs the conservation of both mass 
and angular momentum can be express in the form (e.g. 
iPringldmmh , 



^S=— /, with f = l-J—, 
in \ R 



(13) 



where M is the mass accretion rate, R[ n is the radius of 
the inner edge of the disc taken to be the marginally stable 
orbit, and the (density- weighted) mean kinematic viscosity 
i'(R) is defined by 



7X(R) = f ° ^(R,z)d 

J -2(1 



where 
E(i?) = 



= Q r? 



32 



p{R,z)dz = t^Pozo, 



/ z OCIfj 
p(R,z)dz = — p z 



(14) 



(15) 



(16) 



are the surface density and vertically integrated pressure, 
respectively. 

The final equation comes from energy considerations. It 
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is traditional to assume that the disc is in radiative equilib- 
rium. However, convection, turbulence or other motions can 
also contribute to the transport of energy to the disc sur- 
face. In fact, if the disc has a polytropic vertical structure 
with s — 3, the radiative diffusion law gives us P ra< j pro- 
portional to z. If F = Frad, i.e., if radiative diffusion carries 
the entire energy flux, then JTT} gives us a dissipation rate 
per unit volume independent of z, which is neither consis- 
tent with the a viscosity prescription assumed here ([5]) since 
p = Po(R) [l — (z/zo) 2 ] , nor is in agreemen t with simula - 
tions of radiation-pressure dominated discs ijTurneri 12001 ). 
Therefore, here we use the energy balance equation in the 
form 

J q- z (F rad + F«t m ) dz = J ti(m'ydz, (i7) 

where the energy transport by other motions rather than by 
radiative diffusion is encompassed in F cxtra . Since radiation 
is supposed to carry the entire heat flux at the photosphere 
z = zq, we assume that the extra term integrates to zero so 
that the energy balance can be written as, 

F laA (R,z ) = \{m') 2 (a^j , (18) 

since F ra d (R, z = 0) = by symmetry. 

Here we take /i m = 0.615, kt = 0.33 cm 2 g _1 , and use 
non-dimensional parameters to represent the dependence in 
mass, accretion rate and radius. The radial structure of the 
disc can be represented in terms of E(P), P(R), P{R) and 
H(R), the (density- weighted) scale-height of the disc defined 
by: 

, r° pz 2 dz z i 

H = ~r° —a = T" ( 19 > 
J- 20 P dz 9 

The equations read: 

E = 2.5 x 10 5 m 3/5 a- 4/5 m 1/5 / 3/5 r- 3/5 (l + /3)" 4/5 , (20) 

P = 3.0xl0 22 ma _1 /^ 3/2 , (21) 

(3(1 + /3)~ 3/5 = 3.6 x 10 2 m 4/5 a 1/10 m 1/10 / 4A V" 21/2 °, (22) 

H = 1.7 x 10 3 m 1/5 a- 1/10 m 9/10 / 1/5 r 21/2 °(l + Pf'\ (23) 

where m is the mass in units of Mq, m is the accretion 
rate in units of the Eddington accretion rate assuming an 
efficiency of 0.1, r is in units of the gravitational radius 
and / = 1 — \J r\ n /r, where rin is the dimensionless, spin 
dependent, radius of the marginally stable orbit; H, E and 
P are in cgs units. It should be noted that in the limits 
/?> 1 (P r dominates) and j3 <C 1 (Pg dominates) we recover 
the same dependencies in m, m, a, r and / as in regions 
(a) and (b), respectively, of the Shakura & Sunyaev model: 



(a) p > 1: 

E = 1.9a _1 m~V 3/a /~ 1 , (24) 

P = 2.4 x 10 6 mV^m 1 / 4 / 2 /-" 21 / 8 , (25) 

H = 6.1 x 10 5 mm/. (26) 

(b) /3< 1: 

E = 2.5 x 10 s Q- 4/5 m 3/5 m 1/5 r- 3/5 / 3/s , (27) 



13 n r w in 2 ^,4/5 1/10_1/10 r4/5 -21/20 / 00 \ 

p — 3.6 x 10 m ' a m ' j r ' , (28) 

nr . »„ < n 3 -1/10 • 1/5 9/10 21/20 rl/5 / nr ,\ 

H = 1.7 x 10 a m ' m 1 r ' f ■ (29) 

The expression for the total pressure in both regions is the 
same and is given by (|21|l . The differences in the numerical 
factors are due to the treatment of the vertical structure, to 
the small difference in the values used for kt and /i m , to the 
different normalization of radius, and also because here we 
consider the r<f> component of the stress tensor, T r( f, = fj,r£l' , 
to be — (3/2)ap, instead of the usual Shakura & Sunyaev 
value — ap. 

In order to get a smooth transition between the two re- 
gions, and to allow for a combination of gas and radiation 
pressure throughout the disc, we solve equation (|22[) numer- 
ically. The resultant function, /3(r), is then used in (|23[) and 
(|20l) . in order to get expressions for H(r) and E(r). The ra- 
dial disc structure obtained here is used in the study of the 
warp and eccentricity propagation, which is described in the 
following section. 

One reason that we take some care over the basic disc 
model is that it is important to obtain a smooth density 
profile for the disc in order to calculate the propagation of 
the warp and eccentricity. Also, since H 2 appears in the 
equations for those quantities, the factor of 1/9 in equation 
(|19[) makes a significant difference. 

For comparison, we also include in this section the ex- 
pressions obtained for region (a), if the stress scales with 
the gas pressure instead of total pressure (^ = ap g /fl). The 
semi-thickness of the disc is independent of the viscosity 
prescription, so that (|26[) is still valid. The surface density 
is given by (|27[) while j3 and P are given by: 

P = 3.9 x 10 27 m 13/5 cr 4/5 m 1/5 / 13/5 r-- 18/5 , (30) 

= 1.3 x 10 5 m 8/ V /5 m 1/5 / 8/5 r- 21/1 °. (31) 

In Fig. [2] we show how /3, the surface density, the ver- 
tically integrated pressure and the isothermal sound speed 
(squared), c 2 = po/po = |Sl 2 i/ 2 , vary with radius for a disc 
with a — 0.1, a — 0.5 and m = 10, for two different values 
of the accretion rate. The differences in the disc structure 
when the stress scales with the total pressure or with the 
gas pressure are more evident in the innermost region of the 
disc, and for higher accretion rate, i.e., where radiation pres- 
sure is more important. The exception is the sound speed, 
proportional to the thickness of the disc, since H in the 
radiation-pressure dominated region is independent of the 
viscosity prescription. 

Finally, it should be noted that the model presented in 
this section is Newtonian. The effect of the black hole ap- 
pears only in the choice of the inner radius, which is taken 
to be the marginally stable orbit and is, therefore, depen- 
dent on the rotation of the compact obje c t. Th e relativistic 
correction factors of iNovikov fc Thorng (|l973l ) are not of 
great importance here since the most important relativistic 
effects in the study of warp and eccentricity propagation are 
the apsidal and nodal precession. 
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Figure 2. Log-log plot of the variation of (a) /3 = P r /P g , (b) surface density, (c) vertically integrated pressure and (d) sound speed in 
units of c with radius for rh = 0.1, 0.5, a = 0.1, a = 0.5 and m = 10. The full lines represent the disc model with a combination of gas 
and radiation pressure when the stress scales with total pressure, while the dashed lines correspond to solutions in regions (a) and (b) 
when the stress scales with the gas pressure only. 



3 STATIONARY PROPAGATION OF WARP 
AND ECCENTRICITY 

In the general case, the warp tilt and the eccentricity are not 
only functions of space but also of time, and th eir e volution 
has b een studied by (e.g.) iLubow et al.l (|2002l ) and lOgilviei 
(2001). respectively. However, the characteristic time-scale 
for the precession of a global warp or eccentricity is small 
when compared to the orbital frequency in the binary sys- 
tem, which implies that the frequency of these global de- 
formations is negligible compared to the characteristic fre- 
quencies in the inner part of the disc. Therefore, we choose to 
study the steady shape of a warped or eccentric disc around 
a black hole. 



3.1 Equations 

In this section we present the equations that are used to 
describe the stationary wave-like propagation of warp and 
eccentricity in the disc model presented in the previous sec- 
tion. 

The variation of the tilt of the di sc, W , with radius can 
be obtained from (|Lubow et al.|[2002l l 

PR 3 Q. 2 \ dW ~ 
Q 2 - k 2 + 2ia w Q 2 ) ~d~R 



_d_ 

dR 



This equation describes how propagating bending waves 
communicate the warp through the disc and is valid for 
small-amplitude warps; W describes the amplitude and 
phase of the inclination of the disc. The local wavelength of 
these bending waves is approximately the one given by the 
dispersion relation fl]) for (u),m,n) = (0,1,1), if aw = 0. 
The warp propagation is subject to viscous decay, which is 
described here by a dimensionless viscosity parameter des- 
ignated aw- In the general case of non-isotropic viscosity, 
aw / a as the former is related to the T rz and T^ z compo- 
nen ts of the stress tens or T while the latter is related to T r $ 
fsee lLubow et al1l2002l . and references therein). 

The stationary propagation of a small eccentricity 
thro ugh the disc can be describ ed, in the simplest case, by 
(e.g. iGoodchild fc OrilvidlioOel 



_d_ 

dR 



(7 - ia E ) PR 



, dE 
dR 



+ R 2 ^-E = ER 3 {k 
dR 



D. 2 )E,(33) 



Q 2 )W. (32) 



where E is a (possibly) complex function representing the 
amplitude and phase of the eccentricity at a given radius; 7 is 
the ratio of specific heats. In a strictly isothermal disc 7 = 1 
and the global modes described by this equation have a local 
wavelength which is approximately the one given by the dis- 
persion relation |T} for (u),m,n) = (0, 1,0), if oe = 0. The 
equation describing the eccentricity propagation is based on, 
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and agrees with the local di spersion relati o n of, a 2D disc; 
3D effects are discussed by lOgilvid (|200ll . 120081) . The pa- 
rameter oe in equation (|33|) is essentially a bulk viscosity. 
Effects of shear viscosity are not included as its direct im- 
plementation rnayleadto growing eccentric wav es (viscous 
overstabilitv, iKatol ll978T ) (see also lOgilviel 1200 ll , and refer- 
ences therein). The process of turbulent eccentricity damp- 
ing is poorly known and the simplest way of describing it is 
by using a dimensionless bulk viscosity parameter. 

In order to include relativistic effects in the problem, we 
take = (GM /c 3 )(r 3 ^ 2 + a) and use expressions (|3]) and 
(f3| for the radial and vertical epicyclic frequencies, respec- 
tively. For the surface density £ and vertically integrated 
pressure P, expressions (|20|l and (|21[) are used. 



3.2 Numerical method 



We solve equations (132 p and (|33[) numerically, using a 4th or- 
der Runge-Kutta method. In the case of the warp we use the 
boundary condition that corresponds to zero torque at the 
inner edge: dW/dR(R ln ) — 0. In order for the amplitude of 
the solution to be fixed, we specify the value Wo = W(Rh-i) 
at the inner edge, which corresponds to the (small) inclina- 
tion of the inner edge of the disc with respect to the equator 
of the black hole. For the eccentricity we use similar bound- 
ary conditions: dE/dR(Ri n ) — and Eo = E(Ri n ). To avoid 
the formal singularity of the equations at the marginally 
stable orbit, here we take Ri n = R ms + 5, where 5 <C Ri n 
is an arbitrary value; the solutions for W(R) and E(R) are 
practically independent of the choice of 5. The solutions are 
linear and can be renormalized to obtain any desired Wout 
and i?out , the values of disc tilt and eccentricity at the outer 
radius. This normalisation is particularly meaningful in the 
case of W(R) since W ou t is the (constant) disc tilt at large 
radius, which can be related to observations. 

To solve the warp and eccentricity equations we con- 
sider a black hole with 10 solar masses, i.e. m = 10, and 
spin a = 0.5. The ratio of s pecific heats is given by (e.g. 
iKippenhahn fc Weigertlll996l rl. 



5 + 40/3 + 32/3 2 
3 + 27/3 + 24/3 2 



(34) 



In regions (a) and (b) this expression gives 7 ~ 4/3 and 
7 ~ 5/3, respectively. We regard the dimensionless vis- 
cosity parameter a as con stant throughout the disc as in 
IShakura fc Sunvaevl (|l973T ). Although si mulations suggest 
values of a of the order 10 -2 , according to lKing et al.l (|2007l ) 
observations indicate a typical range of a ~ 0.1 — 0.4. Here 
we choose to fix the viscosity parameter to 0.1. Parameters 
aw, «e and rh can be varied. To mimic the transition be- 
tween different black hole states we fix the viscous damping 
and vary the mass accretion rate. The results and corre- 
sponding discussion are presented in the next section. 



2 Note that in many references /3 is used to represent the ratio of 
gas pressure to total pressure while here /3 = p r /Pg- 



4 RESULTS AND DISCUSSION 
4.1 Undamped propagation 

In Fig. Owe show the variation with radius of the warp tilt 
and eccentricity for increasing accretion rate when aw = 
= tiE. In this case no dissipation is present in the equations 
for global deformations, so they propagate everywhere with 
non-negligible amplitude. Since the imaginary part of the so- 
lutions is zero, the warp comprise s a pure tilt, i.e., the disc 
is not twisted (|Lubow et al.ll2"002h . For small accretion rate 
the warp has an osci l latory shape as previously predicted by 
llvanov fc Illarionovl (|l997h . The wavelength increases with 
radius and, at large R, W tends to a constant value, the 
inclination of the outer disc with respect to the equator of 
the black hole. Note, howeve r, that contrary to the results of 
iBardeen fc Pettersonl (|l975h the inner disc is not necessarily 
aligned with the mid-plane of the com pact object, as pre- 
viously shown by iLubow et al.l (|2002l ). The simulations of 
iFragile et al. I (|2007h also found hints of an oscillatory struc- 
ture characterised by a long-wavelength with no indication 
of the Bardeen-Petterson effect. Eccentricity also has an os- 
cillatory structure but its wavelength is shorter (yet much 
longer than the semi-thickness of the disc) and has a slower 
increase with radius. The shortening of the wavelengths at 
small radius of both the warp and eccentricity implies a 
sharp increase of the gradients dW/dR and dE/dR in the 
inner region of the disc. This is an important feature in 
terms of wave coupling in this region since the growth rate 
for inertial modes that are trapped below the maximum of 
the epicyclic frequency, and interact with global deforma- 
tio ns, is proportional to |dW Vd-R| 2 or \dE/dR\ 2 , as found 
bv lFerreira fc Ogilvid (|200Sl) . 

The basic behaviour of the numerical solutions pre- 
sented in this section can be explained using WKB theory. 
In particular, it is possible to predict how the wavelength 
and amplitude of the several solutions scale with radius. 
A WKB analysis of the secular theories can be performed 
by assuming the deformations W and E to be proportional 
to exp(Jiodi?), where b is a function of R to be deter- 
mined. The dP/dR term in equation (|33|) can be ignored 
in the WKB analysis since relativistic precession dominates 
the propagation of eccentricity. Letting Dw = W(R) and 
De — E(R), we can write equations (|32[) and (|33|l in the 
form 



dR V^dfl) +/il A= °' 



W,E. 



(35) 



The WKB solution is valid where the radial variation of gi 
and hi is slow and is given by, 



A 



(gih, 



v-l/4 



exp 



ly/hi/gidR 



(36) 



In the inviscid case (aw = = oe), gw ~ PR 3 /6r , 
h w ~ T,R 3 D, 2 4ar' 3/2 , g E = 'yPR 3 , and h E « E.R 3 fi 2 6r _1 , 
keeping only the lowest-order terms in the expressions for 
O. 2 - k 2 and ft 2 - Q 2 . Remembering that P = Q 2 H 2 T, (from 
[SJ, we can write 



W oc r 1/8 (EH)~ 1/z exp 



V2, 



24a 

5/2#2 



dR 



£ocr 1/4 (Sif)" 1/ ^exp 



V2, 



jrH- 



■dR 



(37) 



(38) 
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Figure 3. Radial variation of (a) the warp tilt and (b) eccentricity normalised to their values at the outer radius for aw = = c?e 
(i.e., no viscous damping) for (1) rh = 0.2, (2) rh = 0.4, (3) rh = 0.8. The full line represents the real part of the disturbance while the 
imaginary part is represented by the dashed line. A logarithmic scale is used for the x-axis. 



where we have ignored the variation of 7 with radius. If 
the stress scales with the total pressure, these expressions 
indicate that the amplitude of the warp should be propor- 
tional to r~ 5//s in region (a) and to r~ 1//10 / _2//j in region (b); 
the amplitude of the eccentricity is proportional to r-~ 1//2 in 
the radiation-pressure dominated regime and to 7- 1//40 /~ 2/ ' 5 
when gas pressure dominates. If the stress scales with the gas 
pressure only, the amplitudes in region (a) are proportional 
to r 17/40 /" 4/5 and r 11/20 f~ 4/5 in the case of the warp and 
eccentricity, respectively. The results for region (a) can be 
confirmed in Fig. [4] In region (b) r 3> n n , / ~ 1, so that the 



radial variation of both warp and eccentricity amplitudes is 
slow, as seen in Fig. [3] The eccentricity does not tend to 
a constant at large radius because its amplitude is slowly 
increasing when r increases. The warp amplitude decreases 
with r but does not tend to zero as indicated by the WKB 
scaling. The reason for this apparent contradiction is in the 
failure of the WKB approximation at large radii where the 
wavelength tends to infinity. 

The radial variation of the wavelength of the deforma- 
tions can also be understood in light of the WKB theory. 
From equations (|37I38[) , we see that the wavelengths are 
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Figure 4. Solutions for (a) warp and (b) eccentricity propagation in region (a) (/3 2> 1) for aw = = c»e, rn = 10, a = 0.5, rh = 0.1, 
when the stress scales with the (1) total pressure and the (2) gas pressure. The full line represents the real part of the disturbance 
while the imaginary part is represented by the dashed line. The dot-dashed curves show the radial variation of the amplitude of the 
deformations, as predicted by the WKB analysis. 



given by Aw ~ Hnr 5/l4 '/\i r 6a, and Ae ~ 2Hnr 1 ^ 2 ^7 /6, for 
the warp and eccentricity respectively. These dependencies 
are in agreement with the results shown in Fig. [3] where it is 
clear that the wavelengths of the deformations increase with 
radius and with the disc thickness (or equivalently with the 
accretion rate, cf. (|26I29[) ~). The increase in A when rh in- 
creases is evident: in particular, for rh = 0.8 the wavelength 
of the warp becomes so large that it practically loses its os- 
cillatory character, but still maintains the sharp increase in 
dW/dR. Even in this case, the inner d isc is not aligned with 
the eq uator of the black- hole as in the lBardeen fc Pettersor] 
(|l975h picture. Unfortunately the behaviour of the solutions 
very close to the marginally stable orbit cannot be trusted in 
detail because the singularity of the disc model there ought 
to be resolved by the transition to a supersonic inflow. 

4.2 Damped propagation 

In Fig. [5] we show the variation with radius of the warp tilt 
and eccentricity in the case where a significant attenuation 
is present. The values of aw = 0.15 and cve = 0.25 are 
both larger than the viscosity value a used to model the 
disc, and are chosen so that the warp and eccentricity have 
negligible amplitude at the inner region of the disc, for small 
accretion rate (rh — 0.2). For dissipation values smaller than 



these ones, the global deformations propagate everywhere 
with non-negligible amplitude even for small accretion rate, 
although with larger amplitude in the inner region for larger 
/if. 

An important difference between the results presented 
here and the ones in Fig. [3] is the existence of a non-zero 
imaginary part in the solutions. This shows that the reflec- 
tion of the global deformations from the stress-free boundary 
is no longer perfect since they are affected by viscous atten- 
uation and, in consequence, their amplitude is reduced as 
they propagate to smaller radii; a wave with inward group 
velocity, as opposed to a standing wave, is set up. 

It is evident from Fig. [5] that when the accretion rate 
increases the global deformations propagate to the inner re- 
gion of the disc more easily, i.e., the amplitude of warp and 
eccentricity in the inner region increases with rh. This result 
can be explained in light of the WKB analysis introduced in 
4.1. If a small viscosity is present (aw, cie 7^ 0), the WKB 
solution for the warp and eccentricity includes an attenua- 
tion factor due to the presence of a imaginary wavenumber 
fa, 

1 Warp a w / 6a Ecc a E /~6~ , . 
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Figure 5. Radial variation of (a) the warp tilt and (b) eccentricity normalised to their values at the outer radius for aw = 0.15 and 
ag = 0.25 for (1) rh = 0.2, (2) rh = 0.4, (3) rh = 0.8. The full line represents the real part of the disturbance while the imaginary part 
is represented by the dashed line. A logarithmic scale is used for the x-axis. 



If H/r varies slowly with radius, the logarithm of the atten- 
uation factor is 

J \ki\dR~j^ {an = aw, Ob), (40) 

indicating that, for fixed viscosity, the attenuation is smaller 
if the thickness of the disc (or equivalently rh) is larger. 

The results presented in this section indicate that global 
deformations can reach the inner region of accretion discs 
under a wide variety of conditions. In particular, even when 
subject to large viscous attenuation, both warp and eccen- 



tricity can reach the inner disc provided the black hole is 
accreting mass at a large enough rate. 

Some caveats accompany the solutions obtained for 
W{R) and E{R). The equations used assume that these 
quantities are small enough so that non-linear effects can be 
neglected. They also assume that the quantities 1 1 — Q 2 . /Q, 2 1 
and |1 - Q, 2 JQ 2 \ are smaller than, or of the order of, H/r 
which is not true of the inner region of relativistic discs. 
Another important caveat relates to the effects of viscos- 
ity in the propagation of global deformations. Here they are 
parametrised using viscous coefficients aw and oe; this is 
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the simplest way of describing the poorly understood pro- 
cess of turbulent damping in accretion discs. In the partic- 
ular case of the eccentricity equations only a bulk viscosity 
is present so that the complications of viscous overstability 
can be avoided. However, the latter may be present result- 
ing in growth (as opposed to decay) of eccentricity at small 
radii. 



final results. Also, since the terms on the right-hand side of 
equations (|32I33[) dominate the warp and eccentricity prop- 
agation, it is more important to introduce relativistic effects 
in these equations by using expressions ([IJ-Q for the char- 
acteristic frequencies, as they correctly describe apsidal and 
nodal relativistic precession rates, related to Q 2 — Q 2 and 

k 2 -n 2 . 



4.3 Relation to global modes 

The above results are for strictly stationary deformations 
such as may be induced by a misalignment between the ro- 
tational axis of the binary orbit and the spin axis of the 
black hole. We have also computed slowly precessing global 
eccentric modes in binary stars by a method similar to 
iGoodchild fc Qgilviel (|2006h but including the relativistic ex- 
pressions for the characteristic frequencies. When relativis- 
tic precession is included, the global modes have almost the 
same precession and decay rates as in a Newtonian model, 
while the solution in the inner part of the disc resembles the 
stationary solutions described above. 



5 CONCLUSION 

In this paper we studied the propagation of warp and eccen- 
tricity in discs around black holes to determine the condi- 
tions under which these disturbances can propagate to the 
inner regions of accretion discs. High-frequency QPOs have 
previously been identified with inertial oscillations trapped 
in the inner region of discs, and are detected mainly when 
black holes are in the very high state where accretion rate 
is maximum. We find the accretion rate to have a vital role 
in the damped propagation of global deformations. Our re- 
sults suggest that the activation of the inner region, and 
consequent excitation of trapped oscillations by these dis- 
turbances, may be possible only when the accretion rate is 
close to its Eddington value, i.e., when the black hole is in 
the very high state. 

When the propagation of global disturbances (found to 
have an oscillatory structure in the radial direction) is not af- 
fect by viscous damping, the increase in mass accretion rate 
gives rise to a lengthening in their wavelength, in agreement 
with the WKB analysis of the warp and eccentricity equa- 
tions. The most interesting results of our calculations are 
obtained when the more realistic situation of damped prop- 
agation of global disturbances is considered. In this case, the 
increase of accretion rate facilitates the propagation of warp 
and eccentricity, i.e., their amplitudes in the inner region 
increase with m. In particular, when the accretion rate is 
only a small fraction of Eddington and the viscous damping 
is strong enough to completely suppress the propagation of 
global deformations, an increase in rh to the values expected 
in the very high state results in their amplitude in the inner 
region being increased to a non- negligible value. 

A limitation of the disc model used in our calcu- 
lations and introduced in section 2 is the fact that it 
does not include relativistic effects. For example, the law 
for angular momentum conservation in relativistic discs 
l|Novikov fc Thornelll973l ) is slightly different from the one 
used here (cf !13[) . However, relativistic corrections are small 
and their inclusion is not expected to significantly affect the 
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